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Abstract 

A new approach to determine the LHC luminosity is investigated. Instead of employing 
the proton-proton luminosity measurement, we suggest to measure directly the parton- 
parton luminosity. It is shown that the electron and muon pseudorapidity distributions, 
originating from the decay of W^, W~ and bosons produced at 14 TeV pp collisions 
(LHC), constrain the x distributions of sea and valence quarks and antiquarks in the 
range from 3 x 10~^ to ~ 10~^ at a of about 10'' GeV^. Furthermore, it is 
demonstrated that, once the quark and antiquark structure functions are constrained 
from the and production dynamics, other qq related scattering processes at the 
LHC like qq — > W^W~ can be predicted accurately. Thus, the lepton pseudorapidity 
distributions provide the key to a precise parton luminosity monitor at the LHC, with 
accuracies of « ±1% compared to the so far considered goal of ±5%. 
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1 Introduction 



Interpretations of essentially all proposed measurements at the LHC, CERN's 14 TeV proton- 
proton collider project, require a good knowledge of the parton distribution functions at the 
relevant and the collected integrated luminosity. Both omnipurpose experiments, ATLAS 
and CMS [0, consider a luminosity accuracy of ±5% as their goal 0. 

The traditional methods to determine the proton-proton luminosity are size and intensity 
measurements of the beams at the interaction point, as well as event rates of processes with 
previously measured or calculable cross sections like elastic proton-proton scattering 0] and 
QED processes like pp ppe^e~ 0. Unfortunately, clean measurements of the above pro- 
cesses, especially at the high luminosity phase of LHC, are very difficult. Consequently it is 
not obvious that the proton-proton luminosity can be measured with a ±5% accuracy. 

Assuming that the proton-proton luminosity is measured, different experimentally observed 
cross sections are compared with theoretical calculations, using parton distribution functions, 
f(a;, Q^), where x is the fractional parton momentum (x = Pparton/-Ebeam) of the relevant types 
of valence and sea quarks (or antiquarks) and gluons at the considered of the reaction. 
These parton distributions are determined from experimental observables in lepton-hadron 
scattering (DIS processes from fixed target and HERA experiments) and Drell-Yan lepton pair 
production processes at hadron colliders These results, obtained at different Q^, have then 
to be extrapolated to the relevant scale of the studied process. While the x distributions of 
the valence quarks are now quite well constrained, uncertainties for the x distributions of sea 
quarks and antiquarks and gluons remain important. As a result of these structure function 
uncertainties, total cross section predictions of W~^, W~ and Z° boson production at 14 TeV 
pp collisions (LHC), vary currently between 10-20%. Even though the experimental errors are 
expected to decrease considerable during the next years, cross section uncertainties related to 
structure functions will remain important. These uncertainties, combined with the unknown 
contributions from higher order QCD corrections, are usually considered to limit the use of the 
reaction pp —>■ W^{Z^) as an absolute proton-proton luminosity monitor at very high center 
of mass energies . 

The above problems result in luminosity uncertainties, which are larger than the considered 
goal of a ±5% proton-proton luminosity accuracy 0. Consequently, current estimates for the 
achievable accuracies of some measurements at the LHC appear to be somewhat depressing. 
This is especially the case when these uncertainties are compared with the possible small 
statistical errors for many LHC measurements, the current knowledge of quark and lepton 
couplings to the and boson or the already achieved accuracies in high energy e~^e~ 
experiments. 

As a solution for the above problem, we propose a new approach to measure the LHC 
luminosity. This approach is based on: 

1. Experiments at the LHC will study the interactions between fundamental constituents 
of the proton, the quarks and gluons at energies where these partons can be considered 
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as quasi free. Thus, the important quantity is the parton-parton luminosity at different 
values of Xparton [H and not the traditionally considered proton-proton luminosity. 

2. Assuming collisions of essentially free partons, the production of weak bosons, ud — 

— > i'^u, du — i> W~ i~u and uu{dd) ^ Z° — > are in lowest order understood 
to at least a percent level. Cross section uncertainties from higher order QCD corrections 
are certainly larger, but are obviously included in the measured weak boson event rates. 
Similar higher order QCD corrections to other qq scattering processes at different Q^, 
like qq — > iy"'"lV~, can be expected. Thus, assuming that the dependence can in 
principle be calculated, very accurate theoretical predictions for cross section ratios like 
a{pp — > W^W~) I a{pp W^) should be possible. 

3. It is a well known fact that the and production rates at the LHC, including 
their leptonic branching ratios into electrons and muons, are huge and provide relatively 
clean and well measurable events with isolated leptons. With the well known and 
masses, possible x values of quarks and antiquarks are constrained from M^± = sXqXq 
with s = 4:E'^eam- The product XgXg at the LHC {y/s = 14 TeV) is thus fixed to^ 3x 10"^ 
Thus, the rapidity distributions of the weak bosons are directly related to the fractional 
momenta x of the quarks and antiquarks. Consequently, the observable pseudorapidity 
distributions of the charged leptons from the decays of and Z° bosons are also 
related to the x distributions of quarks and antiquarks. The shape and rate of the lepton 
pseudorapidity distributions provide therefore the key to precisely constrain the quark 
and antiquark structure functions and their corresponding luminosities. 

The aim of this paper is to demonstrate the feasibility of this approach and thus improve 
the luminosity measurement at the LHC for quark-antiquark related scattering processes. It 
will be shown that the dynamics of the single weak boson production at the LHC allow to 
constrain the q, q structure functions, the corresponding parton luminosities and therefore also 
the cross sections of other qq related processes. Finally, we suggest that a similar approach 
to gluon related scattering processes might eventually also lead to similar accuracies for the x 
distribution of gluons. 



2 Event rates and the selection of pp W^, W and 

The production of pp — W~^, W~ and Z^ and their identification using the leptonic decays 
have been discussed extensively in the literature 0. In particular, these reactions provide clean 
sources of isolated high pt electrons or muons, and due to their high rate, are often considered 
as a clean and excellent calibration tool at the LHC [|l^ . However, previous studies concluded 
that their use as a luminosity monitor is limited to relative luminosity measurements only [0]. 
The reason for these pessimistic conclusions is based on the predicted cross section variations 



using different sets of structure functions The size of these cross section variations is 



as large as 10-20% as can be seen from table 1. The cross section predictions as well as the 



following simulation results are obtained using the PYTHIA Monte Carlo program ||T2 . 
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Weak boson cross sections (LHC 14 TeV) 




ax BR 


Reaction 


MRS (A) 


CTEQ 2L 


GRV 94 HO 




20.18 nb 


17.32 nb 


21.58 nb 


du^W- ^ tv 


14.24 nb 


12.63 nb 


15.40 nb 


uu{dd) ^Z^ ^ 


3.246 nb 


2.854 nb 


3.456 nb 


qq -> {Z*, 7*) ^ e+e+ {Mfj = 150 - 200 GeV) 


9.71 pb 


8.98 pb 


10.26 pb 


qq^W+W- i+ue-jy 


3.53 pb 


3.30 pb 


3.63 pb 



Table 1: Examples of estimated weak boson production cross sections at the LHC for three 
different sets of structure functions using PDFLIB and PYTHIA programs [O, O . In all cases 
the leptonic branching ratios into electrons and muons are included. 



These cross section variations for single W^, Z^ production are strongly correlated with the 
cross section predictions for other qq related processes. As an example, the corresponding cross 
sections for the reaction qq — > W~^W~ are also given in table 1. Thus, even without looking 
at further details, the uncertainties for multi boson production cross sections at the LHC are 
reduced to about 5-10% if event rates are estimated relative to the production rates of single 
W^, Z^ events. Furthermore such relative measurements reduce also errors from branching 
ratios and detection efficiency uncertainties. 



For the following studies, the MRS (A) structure function set ||TT[ is used as an example 
and reference system. Figure la shows the expected rapidity distribution of and W~ , 
which directly reflect the difference between the x distributions of the u, d valence quarks 
and the sea quark or antiquarks. For small rapidities, corresponding to xi^2 values of 
6 X 10~^, most originate from the annihilation of sea quark-antiquarks and only small 
differences between and W~ are expected. For larger rapidities the originate from the 
annihilation of quarks and antiquarks with very different x values. For example, to produce 
a at a rapidity of about 2.5, one finds the corresponding xi,2 values of the quark and 
antiquark to be xi ^0.1 and X2 ~ 3 x 10~^. As the proton is made of two valence u quarks 
and one valence d quark the production is much more likely than the W~ production at 
large rapidities. 

Figure lb shows the pseudorapidity distributions of the charged leptons originating from 
the decays. Because of the V-A interaction, the differences between the pseudorapidity 
distributions of and i~ especially at large r] values are larger than the ones for the 
and W~ . The reason is that the left handed lepton is emitted preferentially in the 
direction of the incoming quark and the right handed antilepton is emitted opposite to 
the quark direction. Thus the observable charged lepton pseudorapidities reflect not only the 
x distributions of quarks and antiquarks but allow also to some extent a distinction between 
valence and sea quarks at a given x. 

As discussed in the introduction, we want to demonstrate that the dynamics and event 
rates of weak boson production at the LHC accurately constrain the quark and antiquark 
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Figure 1: Expected rapidity (pseudorapidity) distribution of (a) and £^ (b) originating 
from the reaction qq i^u at the LHC {^/s = 14 TeV and the MRS (A) structure 

functions [TTI). The assumed luminosity of 100 pb^^ corresponds to about one day of data 
taking with a luminosity of 10^^ sec~^ cm~'^. 

structure functions and their corresponding luminosity. For this purpose simple event selection 
criteria are used. These criteria closely follow the design characteristics of the proposed CMS 
experiment 0. In detail the following lepton selection criteria are used: 

• Electrons and Muons are required to have pt > 30 GeV within a pseudorapidity of 
\ri\ < 2.40. 

• In order to select only isolated leptons, the transverse energy deposit from other par- 
ticles with Pt > 0.5 GeV and \ri\ < 3, found within a cone R around the lepton 
{R = a/A0^ + Arj'^ < 0.6), should be smaller than 5 GeV. 

^Despite the interest in the very forward region, lepton detection up to much larger |ry| values appears to 
be very difficult. 
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• To reduce possible backgrounds from heavy quark decays and to reject high pt and 
Z° production due to initial state radiation, events with reconstructed jets with Ef > 20 
GeV are removed. The jet momentum vector is reconstructed in a cone R < 0.6 including 
all stable particles, with pt > 0.5 GeV and \r]\ < 3. 

Using these charged lepton selection criteria, pp — > i^u events are required to have 
exactly one isolated charged lepton with 30 < < 50 GeV. The resulting pt spectra of £^ 
and their pseudorapidity distributions are shown in figure 2a and 2b respectively. The used 
kinematic and geometric event selection criteria, result in an event detection efficiency of about 
25% for W+ £+z/, and about 28% for W- tv. 
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Figure 2: The observable (a) charged lepton pt and (b) pseudorapidity rj distributions originat- 
ing from the reaction qq ^ i^u at the LHC (y^ = 14 TeV and the MRS (A) structure 
function including all selection criteria discussed in the text. 

To select events of the type — > Z° — > the presence of a pair of isolated leptons with 
opposite charge {e~^e~ or /i+zi"), with a mi+i- = m^o ± 2 GeV. In addition, the opening angle 
between the two leptons in the plane transverse to the beam has to be larger than 135° and 
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Pt{Z'^) < 20 GeV. As these dilepton events are usually considered to be background free, Z° 
events with large pt can be used to constrain the x distribution of gluons, as will be discussed 
later. 

Events of the type qq — > {Z*, 7*) i^i^ with dilepton masses above 100 GeV have a much 
lower rate. However, these events can be used to study the evolution, up to masses where 
the neutral current sector is well understood (e.g. up to masses of about 200 GeV). At least up 
to these dilepton masses, a measurement of the lepton forward backward charge asymmetry. 



following the method of reference , constrains the ratio of valence and sea u and d quarks 
at different x values. 

Using the above kinematical and geometrical event selection criteria, the efficiency to detect 
both leptons from Z° decays is about 16%, and increases to about 23% for dilepton masses in 
the range between 150-200 GeV. 



3 Sensitivity to the q and q structure functions 

We now study the effects of different structure function parametrisations on the measured 
£^ pseudorapidity distributions {qq — > — > £^z/), and on the reconstructed rapidity 
distribution {qq — > Z° — > 

At the LHC, in contrast to proton-antiproton colliders, the antiquarks have to come from 
the sea. Thus, the pseudorapidity distribution of the positive charged leptons, ud —>■ 
t^v contains the information about the sea d quarks and the valence or sea u quarks. The 
negative charged leptons, du — > W~ i~u carry consequently the information about the 
sea u quarks and the valence or sea d quarks. The rapidity distribution of charged lepton 
pairs, from Z^, {Z*, 7*) —>■ provide the information about the sum of sea u and d quarks 

and the corresponding valence and sea quarks. Consequently, the combination of the different 
observable lepton pseudorapidity distributions should provide some sensitivity to the u, d, it 
and d parton densities over a large x range. 

This sensitivity is first investigated by comparing the weak boson production using two 
quite similar structure function sets, MRS (A) and MRS(H) |[TT|. The main difference between 
these two sets lies in the x parametrisation for the light sea quarks. While the older MRS(H) 
set uses u, d flavor symmetric sea distributions, the MRS (A) set includes a fine tuning of 
the sea quark parton densities with some isospin symmetry breaking, required to describe 
the observation of Drell-Yan asymmetries of Ady = {cTpp — crpn)/{crpp — apn) from the NA 51 



experiment 14 



Figure 3a shows the ratio of a{W^) j a{W~) as a function of the charged lepton pseudora- 
pidity for the two structure function sets. The different parametrisations thus lead, depending 
on the lepton pseudorapidity, to a cross section variation of up to about 10%. The double ratio 
MRS(H)[a(iy+)/a(iy-)]/ MRS(A)[ct(IV+)/ct(W^-)] is shown in figure 3b. The differences of 
about 5-10% between the two sets should be compared with the statistical precision, which 
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is smaller than 1% per bin for an integrated luminosity of only 100 pb~^. Furthermore, both 
sets of structure functions predict almost identical cross sections. The ratio between the 
Z cross sections from the two sets, MRS(H)[(t(Z''^)]/ MRS(A)[(t(Z°)], is also shown in figure 
3b. Combining the obtainable information from W~^, W~ and production, the "fine tuned" 
isospin splitting of u and d sea quarks between MRS (A) and MRS(H) should be detectable 
with good accuracy. 
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Figure 3: a) The detected charged lepton cross section ratio, a{i'^i/)/a{i~h'), originating from 
the reaction qq — > i^u as a function of the lepton pseudorapidity for the MRS(H) and 
MRS (A) structure function parametrisation. b) The relative changes of the charged lepton 
ratios of 3a) between the MRS(H) and MRS (A) parametrisations and also the cross section 
ratio of the Z° production using both parametrisations. 

We have shown that the weak boson rapidity distributions are sensitive to small differences 
between the x distribution of u and d sea quarks and antiquarks. We now go one step further 
and study how well q and q structure function can be constrained from the observable weak 
boson rapidities. For this purpose the different cross sections are studied relative to a 
reference structure function, arbitrarily chosen to be the MRS(A) set. 
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The fraction of weak bosons which are produced from the annihilation of valence quarks and 
low X antiquarks increases strongly with increasing rapidity. The valence quark x distribution 
is already quite well constrained. The main difference between the various structure functions 
comes from the sea q and q parametrisations especially at low x. Thus precise measurements of 
the charged lepton pseudorapidity distributions from decays and the rapidity distribution 
of the Z° events constrain the low x domain of u and d. 

The sensitivity of the measurable lepton rapidity distribution to the different sets of struc- 
ture functions is shown in figure 4 and figure 5. Figure 4 shows the observable ratio of the 
to i~ event rates for three different sets of structure functions and an integrated luminosity 
of 100 pb^^. The difference between the various low x sea quark parametrisations are thus 
reflected in the observable lepton pseudorapidities. Consequently, the shape of the £^ pseudo- 
rapidity distributions provide a strong constraint on the underlying x distribution of quarks 
and antiquarks with x between 3 x 10~^ and ~ 10~^. 
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Figure 4: The ratio of the accepted cross sections a{ud — > — i^u) and (T{dd - 
as a function of the lepton pseudorapidity for four different structure functions 



Figures 5a-c show the ratio of the predicted cross sections from different structure 
functions relative to the reference MRS (A) set. The statistical fluctuations shown in figures 
5a and 5b correspond to the errors from roughly one day of data taking at the initial LHC 
luminosity of 10^^ sec~^cm~^ . The expected event rates are roughly a factor of 10 smaller 
and the errors shown in figure 5c correspond to about 10 days of data taking. 
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Figure 5: Rapidity dependence of the cross section predictions from different sets of structure 
functions relative to the one obtained from the MRS (A) parametrisation; a) for b) for l~ 
and c) for the reconstructed Z'^ . 
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Having demonstrated that the pseudorapidity distributions, originating from weak boson 

decays, are very sensitive to details of the quark and antiquark x distributions one can now 
relate the rate of £^ events in a selected pseudorapidity interval to the quark and antiquark 
luminosity at the given x. Obviously, once the shape of the pseudorapidity distribution is 
accurately known, the event rates need to be measured only for a small pseudorapidity 
interval. For example, counting of events from the process pp — > i^u could be 

restricted to the pseudorapidity range of \r)\ < 0.5. Including all selection criteria one would 
observe roughly 150 000 "clean" luminosity events, corresponding to a statistical error of 0.3%, 
per day at the initial LHC luminosity (~ 100 ph~^ /day). 

Once the quark and antiquark luminosity at 10^ GeV^ and in the x range between 

~ 5 X 10^"^ and 10^^ are determined, accurate cross section predictions of other qq related 
processes should be possible. This is studied for the reaction qq W^W~ . The correlation 
between cross section predictions for single and pair production of weak bosons has been 
pointed out already in section 2 (see table 1). For example, the total cross section predictions 
for the process pp — > between the CTEQ 2L and the MRS (A) parametrisations differ by 
about 15%. However, as we suggest to use the process (y{qq W^) as a reference process, 
one has to relate the cross section of for example (j{qq W^W') to the reference reaction 
(T{qq — > W^). Comparing now the prediction for the relative cross sections between CTEQ 2L 
and the MRS(A) one finds that the difference is reduced to ^ 7.5%. 

As a next step, the parametrisations of the g, q structure functions, especially at low x, 
should be adjusted such that the observed ^^ pseudorapidity distributions are described. As 
the final experimental accuracy for the lepton pseudorapidity distributions will be limited by 
systematics, the limitations of the structure function "fine tuning" are difficult to estimate. It 
is nevertheless worth pointing out that neither the £^ momentum and charge determination 
nor differences between and £^ detection are expected to be problematic. Furthermore, 
backgrounds from different sources and efficiency uncertainties can be controlled by the simul- 
taneous analysis of the and events with isolated electrons and muons. We therefore do 
not expect any principle problem of measuring the shape and the rate of the charged lepton 
pseudorapidity distribution with a ±1% accuracy. Thus even small differences for the sea quark 
parametrisation, like those between MRS (A) and MRS(H), as shown in figure 3a and b, should 
be detectable. One could thus use the difference in cross section for the two sets as a pessimistic 
hmitation of the proposed method. Differences between relative cross section predictions for 
different qq scattering processes and the two parametrisations indicate therefore the size of the 
remaining uncertainties. For example the cross section ratios CF^qq — > W^W~) / a{qq — > W^) 
are 4.74 x 10""^ for MRS (A) and 4.76 x 10"^ for MRS(H). Other qq scattering processes like 
a^qq W^Z^)/a{qq — > W^) show similar stability with predicted ratios of 1.78 x 10"*^ for 
MRS (A) and 1.79 x 10"^ for MRS(H). 

Following the above procedure, i.e. constraining the g, q structure functions and the cor- 
responding parton luminosities, the event rate of weak boson pair production appears to be 
predictable with an accuracy of at least ±1%. 
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4 Outlook and conclusions 



A new approach to the LHC luminosity measurement demonstrates that the x distributions 
of valence and sea quarks and their corresponding parton luminosities can be constrained very 
accurately, using the £^ pseudorapidity distributions from the decay of weak bosons. It is also 
shown that this method leads to very accurate rate predictions of other qq scattering processes. 
For example, the strong correlation between the weak boson pair production and the single 
boson production leads to an estimated experimental luminosity accuracy at the ±1% level. 
This should be compared to the often considered optimistic goal of ±5% accuracy. 

We have not investigated the achievable theoretical accuracies, but believe that many the- 
oretical uncertainties, like the q;s(Q^) uncertainties or still unknown higher order QCD cor- 
rections, contribute in very similar ways to the single and pair production of weak bosons. 
Furthermore, the experimental possibility to measure the x distributions of sea and valence 
quarks and the corresponding luminosities to within ±1% should encourage our theoretical 
colleagues to match this experimental accuracy. 

Finally, we argue that the gluon x distribution and the corresponding gluon luminosity can 
also be constrained in a similar way from accurate measurements of the rapidity distribution 
of gluon dominated scattering processes. In fact, as the g, q luminosity can accurately be mea- 
sured from the weak boson rapidity distribution, the rapidity distribution of gluon dominated 
scattering processes has only to be measured relative to the weak boson rapidity distributions. 
Once the gluon x distribution is known relative to the x distribution of quarks, the weak boson 
rate also provides the luminosity monitor for gluon related signal and background processes. 

The possible experimental accuracy thus depends mainly on how accurate the rapidity and 
distributions of gluon related scattering processes can be measured. A very clean signature 
with a well measurable is much more important than a huge cross section. 

Gluon related scattering processes are gg ^ X and gq{q) — > X. As these processes involve 
jets, measurement problems should be minimized by using processes with small backgrounds 
and well measurable pf Candidates for such processes are high pt events with one or more jets 
and an isolated 7 or a ^ ^'^^~). As the energy and momentum of isolated photons and 
leptons can be measured very accurately, the pt of the jets, assuming transverse momentum 
conservation, can also be determined. Thus, the observables are well measured and should 
provide accurate measurements. 

The production of events with isolated high pt photons or Z*^ are dominated by gq — ^ 'y{Z^)q, 
and qq — > '~f{Z^)g. As shown in table 2, the expected cross sections for these reactions, 
including the branching ratios Z^ —>■ and relatively high pt of 7 and Z^ are still quite 

large. Furthermore, the calculable background corrections from the process qq are expected to 
be small as the cross sections are dominated by the qg scattering process. 

Previous studies of 7-jet final states have shown that jet events with isolated vr'^'s provide 



a considerable background [|T^]. This large background will therefore limit the achievable 
accuracy of such a final state. However, the leptonic Z^ decays provide an excellent signature 
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LHC 14 TeV for MRS (A) and PYTHIA 5.7 


reaction 


cr [pb] 

bOGcV < pt < 100 GeV 


cr [pb] 

lOOGcl' < p, < 200 GeV 


cr [pb] 

Pt > 200 GeV 


qq^Z%^ e+i-)g 


36.4 


6.01 


0.71 


qg^Z\^ £n-)q 


150 


34.8 


4.08 


qq -fg 


717 


74.5 


7.45 


qg 75 


6590 


615 


49.3 



Table 2: PYTHIA cross section estimates for high pt final states of the type )q{9) 
and -fq{g). 

and should allow the selection of essentially background free Z^-jet events. We are not aware 
of any detailed LHC study which demonstrates that this process can indeed be measured with 
accuracies of a few %, but see no obvious reason why the rapidity distribution of the clean 
Z°-jet events can not be measured with an accuracy close to ±1%. 

Unfortunately, the accurate determination due to inherent uncertainties of jet energy 
measurements especially at large rapidities will probably limit the interpretation of the observ- 
able rapidity distribution with respect to the gluon x distribution. Nevertheless, such direct 
measurements of the gluon structure function will provide the highest possible accuracy for 
the X distribution of gluons and might eventually lead to cross section predictions with % 
accuracies for other gluon related scattering processes. 

To summarise, we have shown that the rapidity distributions of and Z^ events at the 
LHC provide directly and accurately the x distributions of quarks and antiquarks. Their rates 
are thus a measure of the corresponding parton luminosities. We have shown that such an 
approach might eventually lead to perhaps 1% accurate cross section predictions of qq related 
scattering processes, like qq — > W^W~ at the LHC. We suggest that the detailed measurement 
of the rapidity distribution of the process qg Zq might provide similar accuracies for the 
gluon structure function and the corresponding gluon luminosity. 



Acknowledgements 

We would like to thank H. Dreiner for many detailed discussions and comments on the 
manuscript. We are also grateful to M. Spira for several detailed discussions and to R. G. 
Roberts for his comments. 



12 



References 



[1] ATLAS Collaboration, W. W. Armstrong et al. , Technical Proposal, CERN Report No. 
CERN/LHCC 94-43, LHCC/P2, 15. Dec. 1994. 

[2] CMS Collaboration, G.L. Bayatian et al. , Technical Proposal, CERN Report No. 
CERN/LHCC 94-38, LHCC/Pl, 15. Dec. 1994. 

[3] See for example reference page 177 and reference page 211. 

[4] TOTEM Collaboration, M. Bozo et al, CERN Report No. CERN/LHCC 93-47 (1993) 
(unpubhshed) and CERN/LHCC 94-39 (1994) (unpubhshed). 

[5] V. M. Budnev et al, Phys. Lett. B 39, 526 (1972) and Nucl. Phys. B63 519 (1973). 

[6] For more details see for example P. N. Harriman et al, Phys. Rev. D 42, 798 (1990) and 
A. D. Martin et al, Phys. Rev. D 50, 6734 (1994). 

[7] See for example SDC Collaboration, E. L. Berger et al. , Technical Design of a Detector 
Report No. SDC-92-201, SSCL-SR-1215 (unpublished) section 3.7.4 and reference [H 
page 177. 

QCD and Collider Physics, R. K. Ellis, W. J. Stirling and B. R. Webber, 1996, Cambridge 
University Press, ISBN 521 58189 3. 

[9] For a recent theoretical analysis and further references therein see for example R.K. Ellis 
et al. , |hep-ph/9704239| FERMILAB-PUB-97/082-T and SHEP-96/37 (unpublished) 

[10] see for example reference ^ page 50 and reference p| page 172. 



[11] For the simulations we have used the PYTHIA 5.7 frame |T^ and the structure functions 
MRS (A), MRS(H), CTEQ2L and the GRV 94 HO as implemented within the PDFLIB. 
For details and further references see H. Plothow-Besch, Comput. Phys. Commun. 75, 
396 (1993); PDFLIB (version 6.06) W5051 CERN Computer library. 

[12] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 

[13] M. Dittmar, Phys. Rev. D 55, 161 (1997). 

[14] NA51 Collaboration, A. Baldit et al, Phys. Lett. B 332, 244 (1994). 

[15] See for example P. Aurenche et al, in Proceedings of the ECFA Large Hadron Collider 
Workshop, Aachen, Germany , 1990, edited by G. Jarlskog and D. Rein (CERN-Report 
No 90-10, Geneva Switzerland 1990), Vol. II. 



13 



